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3 i Abstract. We present a method of generation of the Greenberger-Horne- 

Zeilinger state involving type II and type I parametric downconversion, and 
OO , triggering photodetectors. The state generated by the proposed experimen- 

tal set-up can be reconstructed through multi-mode quantum homodyne 
tomography. The feasibility of the measurement is studied on the basis of 
Monte-Carlo simulations. 



1 Introduction 



A number of proposals for generating the Greenberger-Horne-Zeilinger (GHZ) state 

[|| has been suggested in the literature [gj. Such kind of state is very interesting 

Qv' . as it leads to correlations between three particles in contradiction with the Einstein- 

Podolsky- Rosen idea of "elements of reality" B . In the present contribution we present 
. , a scheme for a complete quantum test of a GHZ state of radiation, not just for a 

'p-ii simple verification of some GHZ correlations, which do not prove that a true GHZ 

state has been produced. In fact, the verification of a state-preparation procedure 

C I needs a complete state-reconstruction technique, whereas correlation measurements ^j 

give identical results for very different states of radiation. In this respect, a crucial 
technique for state-preparation tests is quantum homodyne tomography, in which the 
detrimental effect of non-unity quantum efficiency of detectors is taken into account ab 
initio by the reconstruction algorithms. 

In the following we propose a method for generating a GHZ state through type II 
and type I parametric downconversion, and triggering photodetectors. The proposed 

j^ ■ set-up, although it has low rate of production due to low efficiency for single-photon 

downconversion, however is the only way to generate a "true" GHZ state, without an 
additional vacuum component. The scheme allows a tomographic state-reconstruction, 
whose feasibility here is studied on the basis of Monte-Carlo simulations. 

2 Scheme for the GHZ-state generation 

The scheme for the generation of the GHZ state is sketched in Fig. H. A low-gain 
type-II parametric downconvcrter is pumped by a strong coherent beam to generate 
the state 

|0^(l + 27^)-^/MlO>+7(|l/e,l5o)+e^'^Ml/o,lffe))] , (2.1) 

in the four radiation modes /e,o and 5o,e: where o, e represent ordinary and extraordi- 
nary polarizations, 7 denotes the effective coupling depending on the pump strength 
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r igUre ll Sketch of the experimental set-up for generating the GHZ state. A low-gain |fcy»e-II parametric 
downconverter is pumped by a strong coherent beam and generates the state in Eq. ( |2 . l) ) . Two further 
downconversion processes in type I nonlinear crystals excite modes ao^e, &o,e, Co,e and do,e- A and PBS 
denote a wave plate and a polarizing beam splitter, rieefiectively. The photodetector in the lower part of the 
scheme collapses the field modes in the state of Eq. (12. 51). The photodetector in the upper part checks if the 
second photon in one of the modes /e,o has been spliimptiic first crystal. The resulting state is a mixture of 
a Fock state in modes Co,e and of a GHZ state [see Eq. (p.6[ )], the weights of the two components depending 
on the effective gain of crystal A and on the quantum emciency of detector Dl. 



and the nonlinear susceptibility of the crystal, |0) and |la) represent the vacuum and 
the single-photon Fock state for mode a, respectively. 

The state ^ at the output of the first crystal is then impinged on two type-I nonlinear 
crystal A and B. In this case no classical pump is used and the dyn amics must be 
evaluated without the parametric approximation used to derive Eq. (2.1). For example, 
the unitary evolution describing crystal A is given by 



Ua = exp [xA {alblU + e^^" al&l/o - h.c.)] , 



(2.2) 



XA being proportional to the nonlinear susceptibility of the medium. An analogous 
expression can be written for the crystal B. For simplicity, we assume in the following 
XA — Xb = X- The state at the output of the couple of type-I crystals is given by 

UAUB\o = i^ + ^iY'^'m + 

7[(cosx|l/e) +sinx|lao,l&o))(cosx|l5o) + sinx|lce 14)) + 
e*'^i(cosx|l/o)+e*'^^sinx|lae,16e))(cosx|l5e)+e''^^sinx|lcoldo))]}.(2.3) 

As shown in Fig. yj, a wave plate and a polarizing beam splitter act respectively on 
modes Ce^o and 4,0 according to the unitary transformations 



(4 
{do 



do)/V2 

4)/\/2 



(2.4) 
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When photodetector D2 reveals one photon after the polarizing beam splitter, one is 
guaranteed that two photons have been created in the type-II crystal and that the 
photon impinging on crystal B has been split. This occurs with probability P$ = 
772(7^sin^x)/(l + 27^), ??2 being the quantum efficiency of photodetector D2. The 
corresponding reduced state |<i>) writes 

l*> - -^[cOSx(|l/elCo)-e^(^^+'^-)|l/olCe)) 

+ sinxdlflo Ibo ICo) - e*('^i+'^^+'^«)|lae l&e ICe))] . (2.5) 

Photodetector Dl in the upper part of the scheme monitors the splitting of the photon 
impinging on crystal A. When such detector, characterized by quantum efficiency 771, 
reveals the lack of the "pumping" photon, the resulting output state finally reads as 
the following mixture 

\lCo) P2 = jlr'^'"''^, (2.6) 

I "/ -r^^ 2(1— ?7i cos^ x) 

|lCe) P3^P2 



The overall probability Pg of generating the mixture in Eq. ( |2.6| ) is given by Pg 
P$(l — ?7 i cq s^ y) . One easily recognizes in the ffist component of the mixed state 



in Eq. (2.6) the GHZ state for a suitable arrangement of the phases, namely for 



(fl+ifiA+'PB = TT. 

3 Multi-mode tomographic measurement 

Quantum homodyne tomography is the ffist quantitative technique for measuring the 
matrix elements of the radiation density operator g ^ , which is now used in optical labs 
[[7|. Single-mode homodyne tomography can be generalized to any number of modes. 
However, such a simple generalization needs a separate measurement for each mode, 
which cannot be achieved when modes are not spatially separated. For this reason, in 
Ref. Pi it has been proposed a general method for measuring an arbitrary observable 
of a multi-mode electromagnetic field, using homodyne detection with a single local 
oscillator. Such method is a natural application of a recent group-theoretical approach 
to quantum tomography [pj. In the following we recall the main results, providing 
the rule to evaluate the "unbiased estimator" for a generic (M-t-l)-mode operator. 
The quantum expectation value of the operator can be obtained for any unknown 
state of radiation through an average of such estimator over homodyne outcomes that 
are collected using a single local oscillator by scanning different linear combinations of 
modes on it. The quadrature operator to be measured is given by X{9, V') = [A'' {6, i/j) + 

A{6, ip)]/2 with A{9, ifj) — X]z=o ^~^^'^'^i{S)o-i- The vector u{6) represents a point on the 
Poincare hyper-sphere (for the explicit parametrization see Ref. [g|). By scanning the 
values of V'i G [0, tt] and 9i e [0, 7r/2], all possible linear combinations of modes described 
by annihilation operators ai, with Z = 0, . . . , M, are obtained. The homodyne outcomes 
for XlOjip) can be obtained using a single local oscillator prepared in the multi-mode 
coherent state (^ftol^^) with \ji) = e^^'^ui{9)K/2, where K :§> 1. The expectation value 
for a given operator O is evaluated as follows 

{6) = dp[i'] dp[e] dxp^{x;e,ij)£^[d]{x;e,ij), (3.7) 
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where Pri{x;9,ip) denotes the homodyne probabihty distribution of the quadrature 
X{9,tp) for quantum efRciency 7], and the function £rj[0]{x; 9 , ^p) of x, {9i}, {tfji} has 
the foUowing analytic expression 

M+l ^ + 00 

£^[d]ix;9,^P)^—— die-*+2*^"i*^Tr{d:exp[-2jV^l(^,V)]:}, (3.8) 

jVf! ./n 



with K — 2?7/(2?7 — 1). In Eq. (3.8) we used the notation 



fdf,[i.] = f[f fi, fdp[e]^2''Mll[f dfl, sin2(*^-')+i0,cos0, .(3.9) 



For any given operator O Eq. (^) provides the "unbiased estimator" to be averaged 
over aU homodyne outcomes for the quadrature X{9,ip) of all modes in order to ob- 
tain the ensemble average (O) for any unknown state of radiation. Eq. (3/7) can be 
specialized to the matrix element ({n/}|_R|{?7i;}) of the full joint density matrix. This 
will be obtained by averaging the following estimator B 






/=o 






Ko)t] , 



At/! 
(3.10) 



where fii = max(77i/,ni), i/; = min(m;,n;), and L"{z) denotes the customary general- 
ized Laguerre polynomial of variable z. 



4 Numerical results 



The tomographic measurement on the state in Eq. ( 2^ ) can be suitably performed 
by varying randomly the phases and the polarizations for the couples of modes ao,e, 
bo,e and Co^e and then collecting homodyne outcomes by using three local oscillators. 
Such an experimental arrangement represents an intermediate way of using the multi- 
mode tomographic method of Sect. 3, and the usual metho d bas ed on the product 
of single-mode estimators. By using the estimator in Eq. (3.10) one can measure 
tomographically the expectation value of the projector |0)(0| with 



(|laol6olCo)-|-e'*|lael6elCe)) 



on the state g in Eq. (|2.6|) and compare the result with the theoretical value, namely 



V2 



C{4>)=m\4>) = \vi [1 - cos(0 - V')] 



(4.11) 



We report in Fig. S the resul ts of some Monte-Carlo simulations of the tomographic 
measurement of C{q)) in Eq. ( 4.11 ). 
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Figure 2: Tomographic measurement of C(0) in Eq. ( ft. Ill) for ip = 0. On the loft: r) = 0.85, x = 0.3 ir, 
AT = 6 X 10^. On the right: r] = 0.9, x = OAtt, N = 1.7 X 10'. 



In the simulations the quant um e fficiency of detectors Dl and D2 is rji = rj2 — 30% 
and the phase ip in the state ( |2.6D is (/? = 0. The values of the quantum efficiency 77 of 
homodyne detectors, the coupling x of type-I downconvertors, and the number N of 
simulated homodyne data are reported in the caption of the figures. The results of the 
simulations show that for homodyne detectors with quantum efficiency rj = 85% one 
needs about lO'' data to reconstruct the state with relatively small statistical error. The 
experimental values compare very well with the theoretical ones. The bars represent 
the statistical error, whereas the solid line is the theoretical value of C{4>). In each plot, 
all points are obtained by the same sample of data which causes the evident correlation 
between the statistical deviations. 
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